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ABSTRACT 

The Fornax cluster galaxies NGC 1399 and NGC 1404 are ideal for studying the ef- 
fects of a cluster environment on globular cluster systems. Here we present new optical 
imaging of these two galaxies from both the Hubble Space Telescope's Wide Field and 
Planetary Camera 2 and the Cerro Tololo Inter-American Observatory's 1.5m tele- 
scope. The combination of both data sets provides unique insight on the spatial and 
colour distribution of globular clusters. From B-I colours, we find that both galaxies 
have a broad globular cluster metallicity distribution that is inconsistent with a sin- 
gle population. Two Gaussians provide a reasonable representation of the metallicity 
distribution in each galaxy. The metal-rich subpopulation is more centrally concen- 
trated than the metal-poor one. We show that the radial metallicity gradient can be 
explained by the changing relative mix of the two globular cluster subpopulations. We 
derive globular cluster surface density profiles, and find that they are flatter (i.e. more 
extended) than the underlying starlight. The total number of globular clusters and 
specific frequency are calculated to be N = 5700 ± 500, Sjv = 11-5 ± 1.0 for NGC 
1399 and N = 725 ± 145, Sat = 2.0 ± 0.5 for NGC 1404. Our results are compared to 
the expectations of globular cluster formation scenarios. 

Key words: galaxies: individual: NGC 1399, NGC 1404 - galaxies: interactions - 
galaxies: elliptical - globular clusters: general 



1 INTRODUCTION 

Globular cluster (GC) systems are in some respects (e.g. 
luminosity function), more similar to one another than are 
their host galaxies. This suggests an underlying uniformity 
in the physical conditions under which GC systems and 
galaxies formed. On the other hand, properties such as the 
total number of clusters per unit galaxy luminosity show 
significant variations, and fundamental questions concerning 
the origins of these variations remain unanswered. In par- 
ticular, the roles of tidal encounters, mergers, cooling flows, 
and initial conditions in the formation of GC systems have 
yet to be understood. 

Two galaxies of particular interest in this regard are 



NGC 1399 at the center of the Fornax cluster and a nearby 
cluster member NGC 1404. By studying the GC systems 
of these two galaxies we hope to better understand the en- 
vironmental influences on the formation of GCs and their 
host galaxies. Various estimates of the distance modulus to 
the Fornax cluster have been made in recent years. Here we 
adopt a typical value of m-M = 31.2, which places the clus- 
ter 0.2 mag more distant than Virgo (Jacoby et al. 1992). 
With this distance modulus and no Galactic extinction cor- 
rection, the optical luminosity of NGC 1399 is Mv = -21.74 
(Faber et al. 1989). This cD galaxy is surrounded by hot X- 
ray emitting gas which extends out to at least 38 or 190 kpc 
(Mason & Rosen 1985; Ikebe et al. 1996). The mass-to-light 
ratio increases with galactocentric distance reaching a value 



© 0000 RAS 



2 Duncan A. Forbes, Carl J. Grillmair, Gerard M. Williger, R. A. W. Elson and Jean P. Brodie 



of M/L ~ 70 M0/L0 at ~ 5' (Grillmair et al. 1994a). Well 
within the X-ray envelope, even allowing for projection ef- 
fects, at about 10 (50 kpc) to the south-east of NGC 1399 
lies NGC 1404. It is classified as an El galaxy and has an 
optical luminosity of My = -21.37 (Faber et al. 1989). 

The GC systems of NGC 1399 and NGC 1404 have been 
studied using ground-based imaging by a number of work- 
ers (NGC 1399: Hanes & Harris 1986; Giesler & Forte 1990; 
Wagner et al. 1991; Bridges et al. 1991; Ostrov et al. 1993; 
Kissler-Patig et al. 1997a and NGC 1404: Hanes & Harris 
1986; Richtler et al. 1992). For NGC 1399, Kissler-Patig et 
al. (1997a) gives the total number of GCs to be Ngc = 5940 
± 570. This gives a high specific frequency Sjv = 11 ±4. 
A kinematic study of the GCs around NGC 1399 showed 
that the outer GCs appear to be dynamically related more 
to the cluster of galaxies than they are to NGC 1399 itself 
(Grillmair et al. 1994a). This suggests that accretion of intr- 
acluster GCs (West et al. 1995), tidal stripping (Forte et al. 
1982; Forbes, Brodie & Grillmair 1997) or mergers (Ashman 
& Zepf 1992) may help to explain the high specific frequency 
of NGC 1399. For NGC 1404 the number of GCs is less cer- 
tain. Hanes & Harris (1986) from a photographic study es- 
timated Ngc = 190 ± 80, while more recently Richtler et 
al. (1992) found Ngc = 880 ± 120. The resulting Sjv values 
range from 0.5 ± 0.3 to 2.5 ± 0.3 respectively. The latter 
value is close to the average for ellipticals in the Fornax 
cluster. 

Here we present imaging data from HST's Wide Field 
and Planetary Camera 2 (WFPC2) and the Cerro Tololo 
Inter-American Observatory's (CTIO) 1.5m telescope of 
NGC 1399 and NGC 1404. The WFPC2 data provide ac- 
curate magnitudes and colours with virtually no contami- 
nation from foreground stars or background galaxies. These 
data are complemented by wide field-of-view CCD imaging 
from CTIO. In this paper we will focus on the colour (metal- 
licity) and spatial distribution of GCs. The GC luminosity 
functions in these two galaxies and others in the Fornax clus- 
ter are addressed by Grillmair et al. (1997), and a detailed 
study of NGC 1379, using WFPC2 and CTIO data, is given 
by Elson et al. (1997). 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 HST Imaging 

Details of our HST imaging program are given by Grill- 
mair et al. (1997). Here we use four pointings from that 
programme, namely a central pointing on each of NGC 1399 
and NGC 1404, an outer pointing for NGC 1399 and a back- 
ground field. The outer pointing (F0338) is situated in the 
cD envelope of NGC 1399, on the opposite side of the galaxy 
from NGC 1404. The background field (F0336) is situated 
approximately 1.4 degrees south of NGC 1399 in a blank 
region of sky and serves to measure the surface density of 
background sources. 

Exposure times totaled 1800s in F814W and 5200s in 
F450W for NGC 1399; 1860s and 5000s for NGC 1404 using 
the WFPC2 camera. Three images were taken at one point- 
ing, the other two were offset by 0.5 which corresponds 
roughly to integer pixel shifts in both the PC and WFC 
chips. The images were then aligned and median-combined 



using VISTA software which effectively removed both cos- 
mic rays and hot pixels. The resulting images were found 
to be statistically superior to those reduced using standard 
STSDAS tasks, and enabled us to push our completeness 
limit to about B = 26.5 without incurring large numbers of 
spurious detections. 

Compact sources in each pointing have been detected 
and measured using DAOPHOT II/ALLSTAR software as 
described in Grillmair et al. (1997). For the PC chips the 
galaxy was subtracted off first. Magnitudes have been con- 
verted into standard Johnson-Cousins B and I. No Galactic 
extinction correction has been applied. A colour-magnitude 
diagram for all of the detected sources in each of the four 
pointings is shown in Fig. 1. Completeness tests indicate 
that the 50% incompleteness occurs at B 26.5 for all 
four pointings, with small variations between each point- 
ing. (Note for the PC chips it is about 0.5 mag brighter.) 
In order to avoid any colour bias in our results we have 
excluded all sources fainter than this limit. We have also ex- 
cluded sources brighter than B = 21 as these objects would 
be more luminous than Galactic GCs, and so are probably 
foreground stars. We also applied a colour cut to the data, so 
that only sources with 1.2 < B-I < 2.5 are included. This cut 
is roughly equivalent to -2.5 < [Fe/H] < -1-1.0, assuming the 
Galactic colour-met allicity relation of Couture et al. (1990), 
and covers the range of expected GC metallicities. For NGC 
1399, 95% of all sources lie within this colour selection and 
for NGC 1404 it is 90%. This magnitude and colour selection 
are indicated in Fig. 1 by dashed lines. With these selection 
cuts applied our final sample consists of 572 objects in NGC 
1399 central, 59 in NGC 1399 outer and 208 in NGC 1404. 
The background field has only 14 objects. This background 
field allows us to fairly accurately estimate the proportion 
of bona fide GCs in the three galaxy pointings, i.e. > 98% 
in NGC 1399 central, > 76% in NGC 1399 off-center and > 
93% in NGC 1404. 



2.2 Ground— based Imaging 

Broad-band B and I images of NGC 1399 and NGC 1404 
were taken with the CTIO 1.5m telescope. We used a Tek 
2048 X 2048 array with a pixel scale of 0.44 /pixel, yielding 
a field-of-view 15 on a side. The galaxies were observed 
in 1995 December with typical seeing conditions of ~ 1.5 . 
Data reduction was carried out in the standard way (i.e. 
bias and dark subtraction, fiat-fielding and sky subtraction). 
The total exposure times were 9600s in B and 3660s in I 
for NGC 1399; 9100s in B and 3060s in I for NGC 1404. 
After combining, the images were calibrated using aperture 
photometry from the catalogs of Longo & de Vaucouleurs 
(1983) and de Vaucouleurs & Longo (1988). This procedure 
gave an rms accuracy of ± 0.08 mag for NGC 1399 and ± 
0.05 mag for NGC 1404. No Galactic extinction correction 
was applied. 

Globular clusters were detected automatically using 
DAOPHOT. The detection threshold was set at a conser- 
vative 5(7 per pixel i.e., five times the background noise. 
Selection cuts in the DAOPHOT 'sharpness' and 'round- 
ness' parameters were made to help remove cosmic rays and 
extended objects from the candidate lists. (These parame- 
ters help to remove spikey objects and very non-round ob- 
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jects.) For each detected object we measured a 3 pixel ra- 
dius aperture magnitude and applied an aperture correction 
based on a curve-of-growth type analysis for a dozen isolated 
GCs. The rms from the aperture correction is ~ 0.05 mag. 
Colour-magnitude diagrams for all of the detected sources 
for the two galaxies are shown in Fig. 2. For the CTIO data 
we decided to adopt the same colour selection criteria as 
the HST data, and the same bright magnitude cutoff. For 
the faint magnitude cutoff, we use the 50% incompleteness 
value. Based on an examination of the luminosity function, 
we estimate that for NGC 1399 the 50% value is B ~ 24.3 
and for NGC 1404 it is B ~ 24.2. Thus our CTIO selection 
criteria are 1.2 < B-I < 2.5 and 21 < B < 24.2/24.3, as 
is indicated in Fig. 2 by dashed lines. Even with our selec- 
tion criteria we expect some contamination from foreground 
stars and background galaxies in our samples of 1752 ob- 
jects in NGC 1399 and 858 objects in NGC 1404. This issue 
will be discussed below. We note that although the CTIO 
data only sample the brighter GCs, we do not expect this 
to bias the metallicity distribution. For GCs in the Milky 
Way and M31 (Huchra, Brodie & Kent 1991) and in NGC 
5846 (Forbes, Brodie & Huchra 1996) there appears to be 
no luminosity-metallicity dependence. 



3 RESULTS AND DISCUSSION 
3.1 Spatial Distribution 

Given the magnitude depth and the availability of a back- 
ground field, our HST data are ideal to examine the sur- 
face density distribution of GCs. We first calculate the GC 
density in 5 annuli around NGC 1399. We make use of the 
fact that out to 100 radius the WFPC2 gives almost com- 
plete 180° coverage for one hemisphere of the galaxy (see 
Forbes, Franx & lUingworth 1995 for details). For the NGC 
1399 outer and background fields we simply add up the to- 
tal number of GCs, i.e. 59 and 14 respectively. We have also 
made a small correction (~ 15%) for the expected number 
of undetected GCs at the faint end of the luminosity func- 
tion as determined by Grillmair et al. (1997). Next, we di- 
vide the number of GCs by the appropriate spatial coverage 
to give a surface density. Finally, we subtract off the back- 
ground density of 2.9 objects per square arcmin giving the 
background-corrected surface density of GCs. This is shown 
in Fig. 3. 

Excluding the innermost 2 data points for NGC 1399 
we have fit the data with a function of the form: 
p = Por° 

This fit is shown by a solid line for which po = 126 arcmin"^ 
and a = -1.2 ± 0.2. It is clear from the NGC 1399 outer 
pointing that there are GCs at a projected radius of 8.2 
arcmin from NGC 1399. Extrapolating the profile to 9.7 ar- 
cmin (the projected distance of NGC 1404) suggests that 
there is about one GC per square arcmin associated with 
NGC 1399. This corresponds to half a dozen GCs in the 
WFPC2 field-of-view centered on NGC 1404. The proce- 
dure for estimating the surface density profile of NGC 1404 
is similar to that of NGC 1399, except that we subtract both 
the background density and the expected density of GCs in 
each annulus due to NGC 1399. Again the HST data is fit us- 
ing a powerlaw profile (excluding the innermost data point) 
with po — 36 arcmin"^ and q = -1.3 ± 0.2. 



For both galaxies the GC surface density rises towards 
the galaxy center but flattens off (in log space) in the very in- 
ner regions giving a definable 'core' to the GC system where 
the profile changes slope. This does not appear to be a selec- 
tion effect and has been observed in other large galaxies (e.g. 
Grillmair et al. 1986, 1994b). Forbes et al. (1996) found that 
the 'core radius' of the GC system is loosely correlated with 
the galaxy luminosity. We estimate that the GC systems of 
NGC 1399 and NGC 1404 have core radu of 40" (3.4 kpc) 
and 30 (2.5 kpc) respectively. These values are within the 
scatter of the relation found by Forbes et al. (1996) and are 
similar to the galaxy effective radii. We have measured a GC 
surface density slope of -1.2 ± 0.2 for NGC 1399. There have 
been three previous CCD studies which measured the GC 
slope. Wagner et al. (1991) found that interior to ~ 120 
the GC slope was -1.4 ± 0.11, and this steepened to -1.54 
± 0.15 for larger radii. Kissler-Patig et al. (1997a) derived 
-1.55 ± 0.25 and -1.75 ± 0.3 from two pointings for radii 
beyond ~ 50 . Bridges et al. (1991) measured -1.4 ± 0.2 
from their B band data and -1.5 ± 0.2 from their V band 
data. Our value is somewhat flatter than those measured 
by others. One possible reason for this is that our outer 
pointing is located in the cD envelope of NGC 1399 where 
the GC surface density 'flattens out' (e.g Kissler-Patig et 
al. 1997a). If we exclude the outer point, our fitted slope 
steepens to -1.25. Another possibility is that, even though 
we have not included the inner two data points in the fit, 
we are still fitting into the core-flattened region. The stellar 
profile (log intensity) for NGC 1399 is also shown in Fig. 3. 
We measure a slope of -1.6 ± 0.1 for the galaxy light. This 
compares with the measurements of -1.67 ± 0.12 (Wagner 
et al. 1991) and -1.75 ± 0.1 (Kissler-Patig et al. 1997a). 
There are been some debate as to whether the GC profile 
is flatter than the galaxy profile in NGC 1399. In particu- 
lar. Bridges et al. (1991) claimed that the GC system had 
the same slope as the galaxy. Whereas both Wagner et al. 
(1991) and Kissler-Patig et al. (1997a) said it was some- 
what flatter. Excluding our data, the weighted mean value 
from previous studies is -1.47 ± 0.09 for the GC slope and - 
1.72 ± 0.06 for the galaxy slope. This indicates that the GC 
profile is 0.25 ± 0.11 flatter than the galaxy starlight. Our 
data suggest that the difference is even greater, i.e. 0.4 ± 
0.22. Thus, we conclude that the GC surface density profile 
is fiatter than the galaxy starlight by about 0.3 in the log at 
about the 2 a significance level. For NGC 1404, we measure 
a GC slope of -1.3 ± 0.2 and a galaxy starlight slope of -1.9 
± 0.1. Thus the GC systems in both galaxies have a fiatter, 
more extended distribution that the underlying starlight. 

Using the surface density profile we can estimate the 
total number of GCs and the specific frequency (Sjv; Har- 
ris & van den Bergh 1981) for NGC 1399 and NGC 1404. 
For NGC 1399 GCs are clearly extended beyond 8 arcmin 
in our data. We have decided to use 10 arcmin as the limit- 
ing radius, which is the same limit used by Hanes & Harris 
(1986) in their photographic study. Integrating the power- 
law profile between log r = -0.2 (38 ) and 1 (10 arcmin), 
and adding the number of GCs interior to 38 gives a total 
of 5700 GCs. The uncertainty in deriving the total number 
of GCs is dominated by the choice of limiting radius. We 
have used ±1 arcminute about the limiting radius as our 
error estimate. Thus limiting radii of 9 and 11 arcminutes 
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give N = 5700 ± 500. This is in good agreement with the 
determination of Kissler-Patig et al. (1997a) of N = 5940 ± 
570. For an absolute magnitude of = -21.74, the specific 
frequency Sjv = 11.5 ± 1.0. 

For NGC 1404 the limiting radius is difficult to esti- 
mate. Richtler et al. (1992) determined a background level 
at ~ 200". We find that in the direction of NGC 1399, the 
density of GCs is clearly dominated by NGC 1399 beyond 
about 250 . We have chosen to integrate out to 240 (4 
arcmin) with a reasonable range of 3-5 arcminutes. Thus 
integrating the profile from log r = -0.3 (30 ) to 0.6 (4 ar- 
cmin) and adding the number of inner GCs gives N = 725 
± 145 GCs associated with NGC 1404. This lies between 
the Hanes & Harris (1986) value (190 ± 80) and that of 
Richtler et al. (1992), i.e. 880 ± 120. Compared to these 
ground-based studies, our HST data has the advantages of 
accurate background subtraction (star, galaxies and NGC 
1399 GCs) and nearly 100% complete magnitude coverage. 
For an absolute magnitude of Mv = -21.37, our results in- 
dicate a relatively low specific frequency for NGC 1404 of 
Sjv = 2.0 ± 0.5. Again using the HST data, we calculate the 
'local' Sjv values (i.e. using the number of GCs and the in- 
tegrated magnitude at the de Vaucouleurs effective radius) 
to be ~ 1.5 for NGC 1399 and ~ 0.5 for NGC 1404. 

In Figures 4 and 5 we show the azimuthal distribution 
of GCs, as determined from the HST data, over the range for 
which we have uniform radial and azimuthal coverage. For 
NGC 1399 (Fig. 4) slight enhancements in the GC counts 
can be seen close to the galaxy major axes (position angles ~ 
100° and -80°). There is also a slight deficit around the mi- 
nor ajcis P. A. ~ 10° . Thus the GC system is broadly aligned 
with the stellar isophotes. The direction towards NGC 1404 
(P. A. ~ 150°) is not covered by our HST image. For NGC 
1404, shown in Fig. 5, there is slight deficit in GC counts 
along the galaxy major axis lies at P. A. ~ -20°. There are 
also two enhancements at position angles ~ -50° and 10° 
which do not correspond to either galaxy major or minor 
axes. 



3.2 Colour and Metallicity Distributions 

The GC colour distributions, after colour and magnitude 
selection, for the four HST pointings are shown in Fig. 6. 
Both NGC 1399 and NGC 1404 are dominated by red (B- 
I ~ 2.1) GCs but with a significant tail to the blue (e.g. 
B-1 ~ 1.6). The GCs in the outer pointing of NGC 1399 
are bluer on average than those in the central pointing. The 
median colours are B-I = 1.97 for NGC 1399 in the central 
pointing, 1.70 for the outer NGC 1399 pointing and 1.96 for 
NGC 1404. 

The CTIO data, after colour and magnitude selection, 
are shown in Fig. 7. We do not have a background field 
for the CTIO data. Although we have attempted to remove 
galaxies using the DAOPHOT sharpness and roundness pa- 
rameters, our seeing conditions of ~ 1.5 imply that some 
galaxies will be included in our candidate list. The issue of 
background galaxy contamination in GC colour/metallicity 
distributions is discussed in some detail by Elson et al. 
(1997). Background galaxies in our magnitude range peak at 
B-1 ~ 1.0 and get bluer at fainter magnitudes. The galaxy 
counts are falling fast at our blue cutoff of B-I = 1.2. Given 



the richness of the NGC 1399 GC system, background galax- 
ies are likely to be a small afi'ect. However it may be more of 
a concern for NGC 1404. We can estimate the effects of back- 
ground contamination using the CFRS sample (Lilly et al. 
1995) within the appropriate B magnitude range and scaled 
in area to match that of the CTIO data. The background- 
corrected distributions are shown by a dashed line in Fig. 
7. As we have already removed some galaxies in our initial 
detection process, the true distribution may lie somewhere 
between the two histograms. In either case (uncorrected or 
corrected) the distributions reveal a somewhat different sit- 
uation to the HST data, with both galaxy's GC systems 
peaking in the blue (B-I ~ 1.6) with a significant red tail 
(around B-1 ~ 2.1). The median colours are B-I = 1.84 for 
NGC 1399 and 1.71 for NGC 1404. Thus the CTIO colours 
are on average systematically ~ 0.1 mag bluer than the HST 
data. This result is not formally significant as the error is 
also on the order of 0.1 mag, but differences between the 
shape of the two distributions are clearly seen. As the HST 
data cover only the central regions of each galaxy, this sug- 
gests a radial colour gradient from red dominated objects at 
small galactocentric radii to blue ones further out. In prin- 
ciple, such a gradient could be caused by blue background 
galaxies in the CTIO data. However we give further evidence 
below that this is not the case. 

In Fig. 8 we show the GC colour distributions for NGC 
1399 and NGC 1404 from the HST data. Here the data 
for NGC 1399 includes both the central and outer point- 
ing. The background-corrected histograms are shown by 
dashed lines. Some galaxies, such as M87 (Whitmore et al. 
1995; Elson & Santiago 1996) reveal a clear bimodal GC 
colour distribution. In other cases (i.e. NGC 1374, 1379, 
1387, 1427) the distribution is essentially consistent with 
a single unimodal colour for the GC system (Kissler-Patig 
et al. 1997a,b). The situation for NGC 1399 and NGC 1404 
is not as clear. Nevertheless there is some evidence for multi- 
ple GC populations. Firstly, both galaxies have distributions 
that are significantly broader than that expected from the 
photometric errors (typically a ~ 0.15), i.e. there is a real 
spread in GC colours. Secondly, we have tested the unbinned 
colour data using the KMM statistical test (Ashman, Bird & 
Zepf 1994). This test rejects a single Gaussian fit to the NGC 
1399 and NGC 1404 data with a confidence of over 99%. If 
we represent the data with two Gaussians then KMM deter- 
mines means of B-I = 1.7 and 2.1 for NGC 1399, and B-I = 
1.6 and 2.1 for NGC 1404. Thirdly, as an exercise, we gen- 
erated a colour histogram that is the sum of two Gaussians. 
These Gaussians have mean colours of B-I = 1.6 and 2.1 
with dispersions similar to the photometric error, and with 
total numbers in the ratio of 1:4. They are shown in the 
lower left panel of Fig. 8. In the lower right panel, we show 
the same two Gaussians but with alternating Poisson noise 
included. To the eye, the two noisy Gaussians are qualita- 
tively similar to the HST data, suggesting they could consist 
of two Gaussian-like populations. Fourth, the two peaks in- 
dicated in the HST data are also present in the CTIO data 
(Fig. 7). We conclude that both galaxies do not have a sin- 
gle, uniform GC population but rather show evidence for a 
multimodal GC colour distribution. 

Two other large data sets exist for GCs in NGC 1399. 
They are the CTIO 4m observations by Ostrov et al. (1993) 
using Washington photometry and the Las Campanas 2.5m 
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observations by Kissler-Patig et al. (1997a) in V and I. 
In order to compare all of the different data sets we have 
converted each into metallicity. Broad-band colours can be 
transformed into metallicity using the Galactic relation of 
Couture et al. (1990) given the usual caveat that age effects 
and abundance anomalies may be present. The Washington 
photometry of Ostrov et al. is transformed using the relation 
derived by Geisler & Forte (1990). Washington photometry 
is the most sensitive to metallicity (rms ~ 0.1), then B-I 
colours (rms ~ 0.15) with V-I colours being the least sen- 
sitive (rms ~ 0.25). In Fig. 9 we show the GC metallicity 
distributions from our HST and CTIO data sets, along with 
those of Ostrov et al. and Kissler-Patig et al. None of the 
data sets have been background corrected in this figure (al- 
though for the HST sample the correction is negligible). In 
general, the samples reveal a broad distribution with sev- 
eral minor enhancements, many of which appear common 
to several data sets. From their data, Ostrov et al. claimed 
that NGC 1399 has a trimodal GC metallicity distribution, 
with peaks at [Fe/H] ~ -1.5, -0.8, -0.2. The metal-rich peak 
corresponds to that seen clearly in our HST data. The two 
metal-poor peaks may be real or simply a single popula- 
tion with [Fe/H] ~ -1.0. The mean B-I colours found by 
the KMM test on the HST data correspond to [Fe/H] ~ - 
1.1 and -0.1. The GC metallicity distribution for NGC 1404 
from our HST data is shown in Fig. 10. As indicated by the 
KMM test, a bimodal distribution is a better representation 
of the data than a single Gaussian. The two GC populations 
have means of [Fe/H] ~ -1.5 and -0.1. 

We now return to the CTIO data and the issue of back- 
ground contamination. Earlier we argued that contamina- 
tion, although present, did not change the basic appearance 
of the colour distributions. Further evidence that the two 
peaks in the CTIO data are dominated by bona fide GCs 
comes from examining the spatial distribution of the two 
subpopulations. First, we define the metal-rich and metal- 
poor GC subpopulations in each galaxy as being ± 0.3 dex 
about the mean metallicity of each subpopulation. We then 
calculate the raw surface density (i.e. not corrected for back- 
ground contamination or for the missing faint end of the lu- 
minosity function) for the two subpopulations. These surface 
density profiles are shown in Fig. 11. A background correc- 
tion would tend to lower the metal-poor profile, particu- 
larly at galactocentric radii beyond 4 arcmin in NGC 1404 
(which is excluded from the figure). Although there is con- 
siderable scatter, the surface density of each subpopulation 
declines with distance from its parent galaxy. This indicates 
that each subpopulation is associated with the galaxy and 
is therefore dominated by GCs. Secondly, the fact that the 
HST colours (metallicities) agree with those in the inner an- 
nuli of the CTIO data (see figures 12 and 13) indicates that 
the subpopulations are dominated by GCs. 

Next we investigate the radial variation of GC metallic- 
ity. We have seen that for both NGC 1399 and NGC 1404, 
the metal-rich (red) subpopulation dominates in the HST 
data and in the CTIO data the metal-poor (blue) subpop- 
ulation is dominant. As the HST data probe the inner ~ 
100 region and the CTIO data cover the regions beyond ~ 
30 this would suggest that the metal-rich GCs are more 
centrally concentrated than the metal-poor ones. A simi- 
lar situation is inferred for NGC 4472 (Geisler et al. 1996) 
and NGC 5846 (Forbes, Brodie & Huchra 1997). An un- 



weighted fit to the HST data indeed indicates a metallicity 
gradient in both galaxies. For NGC 1399 we find [Fe/H] = 
-0.31±0.13 log R (arcsec) - 0.046±0.22. This slope is similar 
to that found by Ostrov et al. (1993) i.e. -0.34±0.18. The 
fit for NGC 1404 gives [Fe/H] = -0.36±0.13 log R (arcsec) 
- 0.034±0.22. In both cases the HST fits are consistent with 
the CTIO data. In NGC 4472, Geisler et al. (1996) found 
that the radial GC metallicity gradient was actually due to 
the changing relative mix of the two GC subpopulations. As- 
suming that NGC 1399 and NGC 1404 can be represented 
by metal-rich and metal-poor GC subpopulations we inves- 
tigate the radial gradients further. 

In Figures 12 and 13 we show the metallicity in various 
radial bins for the different GC subpopulations and the mean 
metallicity for the combined population. The range in metal- 
licity included in each data point is ± 0.3 dex. For NGC 
1399 (Fig. 12) the metal-poor and metal-rich GC subpopu- 
lations show no obvious radial metallicity gradient. However 
the mean values for the whole GC system do show a global 
radial gradient. Although not quite as convincing, the same 
trend appears to be true for NGC 1404 (Fig. 13). For both 
galaxies, we conclude that the radial metallicity gradient for 
the overall GC system is consistent with the changing rela- 
tive proportions of the GC subpopulations; the metal-rich 
GCs are centrally concentrated while the metal-poor ones 
are preferentially located in the outer regions. The lack of 
a real abundance gradient places constraints on the role of 
gas dissipation in the formation of GCs. In NGC 4472, the 
metallicity of the galaxy field stars was found to be of the 
same or of slightly higher metallicity than the metal-rich 
GCs (Geisler et al. 1996). For NGC 1399, Brodie & Huchra 
(1991) quote a spectroscopic metallicity of [Fe/H] = -fO.2 
± 0.9. This is similar to the metallicity derived from B-I 
colours of the galaxy (from our CTIO data and from Goud- 
frooij et al. 1994). For NGC 1404, the galaxy B-I colours in- 
dicate [Fe/H] ~ -0.2. Thus for both galaxies, the field stars 
appear to have similar metallicities to those of metal-rich 
GCs (which have [Fe/H] ~ -0.1). 

3.3 Comparison with Globular Cluster Formation 
Scenarios 

3.3.1 The Merger Model 

The merger model of Ashman & Zepf (1992) and Zepf & 
Ashman (1993) make some specific predictions concerning 
the properties of GCs in ellipticals that are the result of a 
gaseous merger. One of their key predictions is that ellipti- 
cals will have multimodal GC metallicity distributions. Fur- 
thermore the metal-rich GCs will be centrally concentrated 
and the more metal-poor ones will be preferentially located 
at large galactocentric radii. This indeed appears to be the 
case for both NGC 1399 and NGC 1404, and is therefore in 
qualitative agreement with their model. For NGC 1404 with 
a low Sjv value (i.e. 2.0 ± 0.5), Ashman & Zepf (1992) expect 
slightly fewer metal-rich (new) than metal-poor (old) GCs. 
Our data are generally consistent with this expectation. For 
NGC 1399 with Sn = 11.5 ± 1.0, they would expect N^/Np 
~ 3-4. In the central regions of NGC 1399 such a ratio may 
hold, but globally N_r/Np < 1. This is of course complicated 
by the difficulty in defining clear metal-rich and metal-poor 
GC subpopulations. However there seems to be a shortfall 



© 0000 RAS, MNRAS 000, 0-| 



6 Duncan A. Forbes, Carl J. Grillmair, Gerard M. Williger, R. A. W. Elson and Jean P. Brodie 



in the number of newly created GCs if gaseous mergers are 
responsible for the metal-rich GCs and the high Sjv value. 
Another expectation from their merger model is that the 
galaxy starlight and GC system should have a similar spa- 
tial distribution in galaxies with high Sat values (i.e. those 
galaxies in which large numbers of new GCs should have 
been created by the merger). Thus we might expect the sur- 
face density of GCs to match that of the galaxy surface 
brightness profile for NGC 1399 but not for NGC 1404. Our 
data indicate that for both galaxies the GC distribution is 
notably flatter than the underlying starlight. We might also 
expect the metal-rich GCs, which should have formed from 
the gas of the progenitor galaxies, to reveal a radial metallic- 
ity gradient as the formation process should be dissipative. 
Figures 12 and 13 suggest that there is little or no metallicity 
gradient for the metal-rich GCs in either galaxy. 

To summarize, the GC systems in NGC 1399 and NGC 
1404 reveal some properties that are consistent with the 
Ashman & Zepf (1992) merger model, but there are also no- 
table conflicts. On an individual basis, some of the disagree- 
ments with the model may be accommodated by slight mod- 
iflcations to the initial model (see Zepf, Ashman & Geisler 
1995). Forbes, Brodie & Grillmair (1997) have recently in- 
vestigated whether a sample of elliptical galaxies meets the 
general predictions of the Ashman & Zepf (1992) model, and 
concluded that mergers were unlikely to account for the GC 
systems in large elliptical galaxies. They favored a multi- 
phase collapse scenario. Next we compare the data for NGC 
1399 and NGC 1404 with this scenario. 

3.3.2 The Multiphase Collapse Model 

Forbes, Brodie & Grillmair (1997) have proposed that GCs 
in ellipticals form in situ during a multiphase collapse. In 
the flrst pre-galactic phase, metal-poor GCs are formed 
and later on, in the second galactic phase, the metal-rich 
GCs form. This leads to a bimodal metallicity distribu- 
tion with the metal-rich GCs more centrally concentrated 
than the metal-poor ones. In massive ellipticals the collapse 
could be largely dissipationless and so no strong radial abun- 
dance gradients are expected. The radial distribution of (the 
metal-poor) GCs will be more extended than the galaxy 
itself, if as expected, they form before the galactic stars. 
The overall metallicity and spatial properties of GCs derived 
from our data are consistent with this picture. 

For cD galaxies, such as NGC 1399, Forbes, Brodie & 
Grillmair (1997) suggested an additional source of GCs from 
tidal stripping of nearby galaxies. They suggested that some 
of GC subpopulation identified by Ostrov et al. (1993) at 
around [Fe/H] ~ -0.8 were acquired from NGC 1404 and 
possibly other Fornax galaxies. Recently, Jone et al. (1997) 
have interpreted the asymmetric X-ray emission around 
NGC 1404 as evidence for tidal interaction with NGC 1399. 
The idea of tidal stripping of GCs is given some support by 
our new data. The outer regions of NGC 1404 are dominated 
by metal-poor GCs, and it is these GCs that may be prefer- 
entially stripped and captured by NGC 1399. Thus some of 
the GCs in NGC 1399 with this intermediate metallicity may 
have originated from NGC 1404. We confirmed that NGC 
1404 has a remarkably low specific frequency value (Sat = 2.0 
± 0.5) for a cluster elliptical, suggesting that some GCs may 
have been lost. If NGC 1404 originally had Sat = 5, i.e. typ- 



ical of normal ellipticals outside of the Fornax cluster then 
it has lost ~ 1000 GCs, which may have been tidally cap- 
tured by NGC 1399. If we assume that the GC system and 
galaxy starlight originally (i.e. before tidal stripping) con- 
tinued beyond 4 arcminutes with the same slope as found 
in section 3.1, then we can calculate the Sjv value of the 
stripped GCs. We find that the stripped GCs are unlikely 
reach values of S]v ~ 11.5 in sufficient numbers to have any 
appreciable effect on the specific frequency of NGC 1399. We 
expected the GC surface density proflle to match that of the 
underlying starlight in the inner regions of NGC 1404, but 
it appears to be somewhat flatter than the stellar profile (at 
least for galactocentric radii of less than ~ 150 ). This may 
just be indicating that tidal stripping has not yet affected 
the galaxy inner regions. The above circumstantial evidence 
suggests that NGC 1399 has acquired GCs from NGC 1404, 
although this is unlikely to explain the high specific fre- 
quency of NGC 1399. Large samples of GC metallicity and 
kinematic information from spectra are probably required 
before the multiphase collapse and tidal stripping ideas can 
be fully tested. 



4 CONCLUSIONS 

We present new optical data of the globular cluster (GC) 
systems in NGC 1399 and NGC 1404. Our images were taken 
in B and I filters using the Hubble Space Telescope (HST) 
and the 1.5m telescope at Cerro Tololo Inter-American Ob- 
servatory (CTIO). This provides us with high spatial reso- 
lution in the central regions and wide area coverage respec- 
tively. From the HST data we have detected over 500 GCs 
in NGC 1399 and over 200 in NGC 1404, covering ~ 90% of 
the GC luminosity function. With the aid of an HST back- 
ground field of the same exposure time, we estimate that 
the contamination levels in the HST samples are only a few 
percent. 

Both galaxies reveal a broad GC colour/metallicity dis- 
tribution. These distributions are inconsistent with a single 
Gaussian. For NGC 1399, two Gaussians with mean metal- 
licities of [Fe/H] ~ -1.1 and -0.1 provide a reasonable rep- 
resentation of the data, although a multimodal distribution 
with several different GC subpopulations remains a possibil- 
ity. In the case of NGC 1404 we tentatively claim evidence 
for a bimodal GC distribution with the metal-poor popu- 
lation around [Fe/H] ~ -1.5 and the metal-rich population 
at [Fe/H] ~ -0.1. In both galaxies the different metallic- 
ity peaks are consistent with distinct GC subpopulations in 
which the metal-rich GCs are more centrally concentrated 
than the metal-poor ones. Globally, the GC systems reveal 
metallicity gradients. However these gradients are consistent 
with a changing relative mix of the two GC subpopulations 
with galactocentric radius. 

We derive GC surface density profiles, and show them to 
be flatter (i.e. more extended) than the underlying starlight. 
The total number of GCs and speciflc frequency are calcu- 
lated to be N = 5700 ± 500, Siv = 11.5 ± 1.0 for NGC 1399 
and N = 725 ± 145, Siv = 2.0 ± 0.5 for NGC 1404. Finally 
we discuss our results in the context of two GC formation 
scenarios. The GC data on NGC 1399 and NGC 1404 are 
generally more consistent with a multiphase collapse (e.g. 
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Forbes, Brodie & Grillmair 1997) than a merger (e.g. Ash- 
man & Zepf 1992) origin. 
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Figure 1. Colour— magnitude diagrams for the compact objects in four HST pointings: central NGC 1399 (circles), central NGC 1404 
(squares), outer NGC 1399 (circles) and a background field (triangles). The dashed box indicates the colour and magnitude selection 
criteria (see text for details). 
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Figure 2. Colour— magnitude diagrams for the compact objects in two CTIO pointings centered on NGC 1399 (circles) and NGC 1404 
(squares). The dashed box indicates the colour and magnitude selection criteria (see text for details). 
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Figure 3. Surface density profiles for the globular cluster systems in NGC 1399 (circles) and NGC 1404 (squares) from HST data. The 
surface densities have had a correction for background contamination and for the missing faintest globular clusters. The globular cluster 
systems of both galaxies reveal a 'core region'. Powerlaw fits to the data beyond the core region are shown by solid lines. The surface 
density of compact background objects in HST images is shown by a long— dashed line. A typical error bar is shown on the left. The 
underlying galaxy starlight profile, arbitrarily normalized in the vertical direction, is shown by a short-dashed line. For both galaxies 
the starlight is more concentrated than the globular cluster system. 
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Figure 4. Histogram of globular cluster position angles within 100 radius of NGC 1399, from HST data. The galaxy major axes are 
P. A. ~ 100° and -80°. Globular clusters appear to be concentrated close to the galaxy major axis and show a deficit along the minor 
axis. 
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Figure 5. Histogram of globular cluster position angle within 100 of NGC 1404, from HST data. The galaxy major axis is P.A. 
-20°. Globular clusters appear to show a slight deficit near the galaxy major axis. 
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Figure 6. Histograms of globular cluster B— I colours for the four HST pointings. In the central pointings, both galaxies are dominated 
by red globular clusters with an extended blue tail. 
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Figure 7. Histograms of globular cluster B— I colours for the two CTIO pointings. The dashed line shows the distribution after correction 
for background contamination. In these wide area pointings, both galaxies are dominated by blue globular clusters with an extended red 
tail. 
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Figure 8. Histograms of the globular cluster B-I colours in NGC 1399 and NGC 1404 from HST data (upper panels). For NGC 1399 the 
central and outer pointings have been co-added. The dashed line shows the distribution after correction for background contamination. 
The lower left panel shows the sum of two Gaussians with peaks at B-I = 1.6 and 2.1. The lower right panel shows the same two 
Gaussians with Poisson noise included. The metallicity distribution of both galaxies can be plausibly represented by two Gaussians. 
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Figure 9. The metallicity distribution for globular clusters in NGC 1399 from four different samples. Data from this paper are shown 
in the left hand side panels. The upper right panel shows the data of Ostrov et al. (1993), and the lower right panel shows the data 
of Kissler-Patig et al. (1997a). None of the samples shown here have been background corrected (although for the HST sample this is 
negligible). The Ostrov et al. data have the highest metallicity sensitivity. The HST distribution is inconsistent with a single Gaussian 
and is better described as multimodal. Common peaks (particularly at [Fe/H] ~ -0.2) can be seen in each data set. 
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Figure 10. The metallicity distribution for globular clusters in NGC 1404 from HST data. A KMM test indicates that the distribution 
is inconsistent with a single Gaussian and that a bimodal distribution with mean metallicities of [Fe/H] ~ -1.5 and -0.1 provides a better 
fit. 
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Figure 11. Surface density profiles for the metal— rich and metal-poor globular cluster subpopulations in NGC 1399 (circles) and 
NGC 1404 (squares) from CTIO data. No correction has been applied for background contamination or missing globular clusters at the 
faint end of the luminosity function. Typical data error bars are shown on the left. In both galaxies the two subpopulations decline 
with galactocentric radius indicating that they are dominated by globular clusters. The metal— rich subpopulation is more centrally 
concentrated. 
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Figure 12. Radial variation of metallicity for the globular cluster subpopulations in NGC 1399. The filled circles show the HST data 
and the open circles the CTIO data. The metal-rich, metal— poor and mean for the whole globular cluster system are shown. The mean 
metallicity gradient is consistent with the changing relative proportions of the subpopulations with galactocentric radius. 
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Figure 13. Radial variation of metallicity for the globular cluster subpopulations in NGC 1404. The filled circles show the HST data 
and the open circles the CTIO data. The metal-rich, intermediate, metal-poor and mean for the whole globular cluster system are shown. 
The mean metallicity gradient is consistent with the changing relative proportions of the subpopulations with galactocentric radius. 
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